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ABSTRACT 

The distribution of scale sizes of the diffraction pattern produced by ionospheric 
irregularities responsible for satellite scintillations has been obtained. They 
have been derived from measurements of the time autocorrelation function of 
the fluctuations of the signals received from synchronous satellites and from 
drift velocities obtained by cross correlation methods usingtwo spaced antennas. 
The sizes have been found to be of the order of 200 and not larger than 400 meters. 

The scale sizes of the irregularities as well as the relative amplitude of the 
fluctuations, for the case of lower satellites, has been deduced from the meas- 
urements made with synchronous satellites using a full wave theory. It is found 
that, as the satellite gets closer to the layer of the irregularities, the scale size 
of the diffraction pattern on the ground gets larger. When the distance between 
the satellite and the layer becomes smaller than a Fresnel distance, the r.m.s. 
of the amplitude scintillations is largely reduced. The preferred heights of the 
layer of irregularities has been determined from radar measurements at 
Jicamarca. 

It is recommended that the minimum distance between antennas for a space di- 
versity reception system should be 300 meters for coverage of high altitude 
satellites and 1000 meters for coverage of all satellites. The maximum distance 
is limited by economical considerations and possible problems of real state. 
There is also a technical problem associated with very large spacing and the 
phase (or time delay) of the information received at both antennas, when the 
bandwidth of the modulation is large. This can be solved by a scheme which 
makes use of the knowledge of the position of the satellite to introduce the proper 
delay. 
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A PROGRESS REPORT ON SCINTILLATION OBSERVATIONS 
AT ANCON AND JICAMARCA OBSERVATORIES 


INTRODUCTION 

The purpose of this report is to present a preliminary analysis of the data 
collected in connection with NASA contract work with ESSA (ESSA contract No. 
E22-9-70(N)). Under this contract a study of equatorial ionospheric irregulari- 
ties using satellite transmissions is being performed. 


To carry out this study, the telemetry system of the Ancon Satellite Tracking 
Station is being used to receive 136.4 MHz transmission from the geostationary 
satellite ATS-3 on two spaced antennas. Whenever it is possible observations 
with the Jicamarca radar are also conducted coincident with those at Ancon. 

The primary objective of this study is to relate these measurements to space 
diversity systems that will permit one to reduce the degradation of satellite 
transmissions during the presence of ionospheric irregularities. The data- 
collected will also be used to improve our understanding of the physical causes 
of these irregularities. 

The first part of this report includes an analysis of propagation through a layer 
of irregularities . The purpose of this analysis is to relate the mentioned meas- 
urements with geostationary satellites with those effects that are found on 
lower satellite transmissions. 

In the second part a review of space diversity techniques is made. Finally, the 
third part includes a description of the data collection work and data processing; 
also, use is made of derivations of part one to analyze the information obtained. 


PROPAGATION THROUGH A LAYER OF 
IONOSPHERIC IRREGULARITIES 

In space diversity systems the separation between antennas is chosen to mini- 
mize the probability of simultaneous fading of the received signals. This is ob- 
tained when the separation exceeds the correlation distance of fading on the 
ground. 

For fadings of satellite signals caused by ionospheric irregularities, the corre- 
lation distances projected on the ground depend on both the altitude 
of the satellite and of the irregularities. 
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Jicamarca observations show that ionospheric irregularities have a layered 
structure. Thus a treatment of propagation of radio waves through layers of ir- 
regularities is required to relate correlation distances observed with a geo- 
stationary satellite to those distances observed by lower satellites. Although it 
is not uncommon to find multiple layers of irregularities, for the purpose of 
this report only one layer is considered as a first approximation to the problem, 
and a monographic study of what is known on propagation through a region of 
irregularities is made. In a future work propagation through multiple layers 
will be studied. 

The problem of propagation through a region containing refractive index ir- 
regularities has been treated by Chernov (1961), Tatarski (1961), Yeh (1962), 
Budden (1965) and others. They arrived at equivalent solutions following differ- 
ent approaches. 

In this monograph the problem of propagation through a layer of irregularities 
will be considered as a single scattering problem. The received signals will be 
considered as the sum of the direct signal, E d , plus the signal scattered by the 
irregularities, E s . Multiple scattering is neglected. The geometry of the prob- 
lem is illustrated in Figure 1. Spherical waves are emitted by a source at P 
and received at Q after travelling through a layer of irregularities. The thick- 
ness of the layer is W and is of infinite extent in the x and y directions. 

The direct signal at the receiver is: 


E - e j " k R ] (1) 

d R 


where 


R = PQ 

and the scattered signal can be computed as follows: From Stratton (1941) the 
far-zone fields produced by an electron with acceleration x is given by 


E * = 


e \ sin 6 .. 
- — ) x 

4 7 T 


(2) 


where 


yu 0 is permeability of free space 
e is the charge of an electron 
r, 6 and <£> are defined in Figure 2. 
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Figure 1. Geometry of scattering from a slab of ionospheric irregularities, 
(a) Transmitter at£m from Z axis, (b) Receiver at ^m from the Z axis. 
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Figure 2. Geometrical parameters of scattering produced by nonrelativistic electron. 

For non-relativistic velocities the acceleration of an electron due to an incident 
wave is mainly due to the electric field of the wave, that is 



Combining (2) and (3) we obtain for a single electron 


(3) 





where 



e 47r m 


Let the irregularities of the layer be characterized by a random variable AN 
(x, y, z) electrons/m 3 with average APf - 0. Then the field scattered by AN dv 
electrons is 
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where 


r s is defined in Figure 1. 


and 


E. 

inc . 



(6) 


If the typical scale size of the irregularities l » then most of the field scat- 
tered will be confined to an angle p = k/l measured from the direction of 
propagation. In the ionosphere i is in the order of 300 m and for a 136 MHz 
satellite signal, A. = 2.2 m, then p will be of the order of 0.4°. Therefore, the 
angle 9 will be 90° + P £ 90° for most of the scattered field and 


**e F'O j I- * 
d E (Q) Z-L-l e 3 


k ( + r s )J 


s 


AN d v 


(?) 


r r. 

S 1 


in integral form 


E 

S 


(Q) = r e E 


o 



k ( r i + r s )] 


d x d y d z 


( 8 ) 


Since the region of interest is confined to a small angle from the z axis in 
Figure 1, the usual approximation can be made.. For the exponent part 


r. 

1 




x 2 + y 2 
2 (Ri +~z) 


r 

s 


% (R s - z) + 


(x - ^) 2 + y 2 
2 (R £ + z) 


where £ is defined in Figure la, and for the denominator 

r. r £ (R. + z) (R - z) 

1 S v 1 ' ' S 



then (8) becomes 


E = r E 

s e 


r an 

d J t 


exp - j k 


y* + x - 


+ z )\ 2 
R; +R / 


2 b 


d v 


( 9 ) 


where 


b = 


( R i + z) ( R s - z) 


R. + R 

1 S 


dv“dxdydz 

Normalizing variables and lengths 


x' = k X 

>> 

& 

ii 

z' = k z 

b' 

“kb 

1! 

3 

R' = k R 

R' = k R 
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y' 2 + \ x 
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( 10 ) 


If the irregularities are assumed to be weak, then the scattered field is only a 
small fraction of the total field and the fluctuations are mainly due to the com- 
ponent of E s in phase with E d , Bowhill (1961). Let 


A = 


R (E ) 

e v s 7 


be the amplitude fluctuation coefficient, then 


A CO) 


, r e Tan 

k 2 J b' 


r 

r ? (R* + z ')i 


Jy ' 2 + 

x' - - — 


l 

R. + R 

l— 1 S -J 

( 

l 2 b ' J 


( 11 ) 
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This expression is similar to the one obtained by Yeh (1962), therefore, the 
results he obtained will be appropriated to obtain the mean signal value and 
correlation functions of the amplitude fluctuations. If the case of Figure 2b is 
considered it is found that 


A (Q) % _ e 
k 2 J 


r J 

AN 1 

<< 

to 

+ 

L, _ e ' 

L + r : j 

■i 

_i 

~ C ° S 

L 2 b' J 


d v 


( 12 ) 


where if z is replaced by -z and y by -y it is found that the resulting expression 
is equal to the one that would be obtained if the source were placed at Q and the 
receiver at P. This could have been inferred from the reciprocity principle. To 
compute the mean square value of the amplitude fluctuations £ can be taken 
equal zero to Figure la and 



P 12 

b i b ; 


,'2 


+ X 


> 2 \ 


r> 2 


+ X 
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2 b; 


co s 


2 b' 


dv ; d v 2 


(13) 


where the bar indicates spatial average 


P 12 = 


ANj A N 2 
AN 2 


In this expression spatial average have been taken as equivalent to ensemble 
average. Equation (13) can be transformed into 

A 2 = 2^ AN 2 (I x + I 2 ) 

where 


I. 


■\ / 


P 


r* 2 


12 


b ; 


b 2 
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+ y? 

b i 
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+ y 
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d v; a v ' 
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12 


4 b' 
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assuming a Gaussian correlation function 


Pl2 = 


A Nj A N 2 
AN 2 


= exp - 


/ x2 y 2 

\ X y 



( 15 ) 


where \ , l y and l z are characteristic lengths, or correlation distances, in 
the x, y and z directions. 

Using results obtained by Yeh (1962) 


I, = W' TT S/ 2 l 

1 z 


(16) 


W' TT 'C' 


■R: + W„ 


R: - W /2 



Cd + D?) 1/2 (1 + D 2 ) + (1 - D D )] 

- l LL dy' (17) 

2 (1 + D 2 ) (1 + Bj) 


where 


y' = (R j + + (R * * 

2 


D x = 


D = 
y 


4 y' (R/ + R' _ y') 

■ 

4 y.' (R! + R' s - y') 


*; a 

The D's are equivalent to the wave parameter defined by Chernow (1961). If the 
layer is considered to be thin the D's will not change appreciable between the 
limits of integral (17) and mean values can be assumed and taken outside the 
integral. Let the mean values be 


D = 

X 


4 R. R 


k V (R. + R ) 
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D 4 R i R s 

y " k (R l + R.) 


then (17) becomes 


r 2 % - w' 77 5/2 


U f D 2 ) (1 + a 2 Djft + ~ aP x ) 1/2 

2 (1 + D 2 ) (1 + a 2 D 2 ) 


where 


a= F 


Finally using (16) and (18) in (13) 


A2 = 2^-— j AN 2 k 2 W 77 5/2 < 1 


(1 + D 2 ) (1 + a 2 ~5I)'+ (1 - ttP,D y > 

2 (1 +D 2 ) (1 + a 2 D 2 ) 


Now, if two receiving antennas at Q, (^/ 9 > °> R s ) and ^2 R / ]„^ re C ° n 

sidered, the correlation function A(Q t ) A (Q 2 ) can be obtained using ( ) 


A (Qj) A (Q 2 ) = I — j AN' jj ^ 


ll\ A N 2 P ] 2 cos ( $ i) ’ COS d v i d v 2 


where 


$ = cos 


• „ , ? <*; + z :>T 

y - + *■ - 2 (r; 7 io J 


m = 1, 2 
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Let 


A (Q t ) A (Q 2 ) 
A (Q) 2 


be the normalized transverse correlation function of A. The same procedure 
used to find A 2 can be used (Yeh, 1962) to find 


P A = 


* 3 + *4 

I 1+ I 2 


(21) 


where 


X 3 = 


/ 


12 


4 b y b' COS d d V 2 


1 2 


X 4 = 


P. 


12 


4W~W C ° S ($ i + V d d v 2 
12 


The solutions of I 3 and I 4 are: 


where 



0 = 


e 


k < s ; + 
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Expressions (22) and (23) show that in general the correlation function, p A , of 
the amplitude fluctuations do not have the same shape as the correlation func- 
tion, P 12 , of the ionospheric irregularities. Two asymptotic cases, D x » 1 and 
D « 1, can be used to illustrate this point. For d » 1 it is obtained 

X 7 x 



AN 2 k 2 W tt 5/2 


(24) 


P A ~ exp 


£ 2 

(M-e x ) 2 


(25) 


where 


For D x « 1 


M = (R. + R s )/R. 


A 2 Z ~ ( -Jl\ AN 2 k 2 Wt 7 s/2 l D 2 
4 \k 2 / 


(26) 


Pa 


1 - 4 




(T 2 

(M^) 2 . 


(27) 


Expressions (25) and (27) are approximately Gaussian and have characteristics 
lengths of the order of Ml x . Intermediate cases produce characteristic lengths 
which do not depart much from the two extreme cases (Bowhill, 1961). This 
permits us to use equations (25) and (27) as a first approximation to the problem 
of finding characteristic_lengths of the ionospheric irregularities from satellite 
transmission, although do not depart much from unity. 

At the magnetic equator ionospheric irregularities that cause amplitude scintil- 
lations of satellite signals are elongated along the earth's magnetic field lines 
with elongation ratios of 30 to 60. The small dimension of the diffraction pattern 
is of the order of 225 m and the irregularities appear at an altitude of the order 
of 300 km (Koster, 1967). If the longest direction is taken along the y axis of 
Figure 1 then an elongation ratio of 30 gives a ~ 10 -3 . For this value of a 
equations (19) and (21) show that A 2 and p K are approximately independent of 
a . This means that scintillations at the equator are mainly due to the dimen- 
sion of the irregularities perpendicular to the earth's magnetic field. This 
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causes the time variations of amplitude fluctuations due to a drif ting layer of 
irregularities to be mainly due to the EW component of the drift. 

Two cases will now be considered, the case of transmissions from a synchronous 
satellite and the case of transmissions from a lower satellite. The parameters 
listed in the first column of Table I will be used. 

For a synchronous satellite the transmitted radio waves are practically plane 
when they are incident on the layer of irregularities . The projected correlation 
distances on the ground are approximately equal to those on the ionosphere. 

This is expressed by the fact that M ~ 1. Then in this case the value of p x can 
be obtained from the correlation function p a . 

The fact that the ionospheric irregularities are drifting will permit one to 
measure p A {£) from the time auto- and cross -cor relation functions, p k (r) and 
P A1 A2 ( T )» This can be shown as follows: Assuming that the drifting irregulari- 
ties do not change appreciable during a time r Q = /v D , where v D is the com- 
ponent of the drift velocity in the EW direction, then the field measured at the 
point Qj (0, 0, R s ) at t = t Q will be the same as the field measured at Q 2 (<f , 0, R s ) 
at time t = t Q + £/v D . These arguments show that for the case of a synchronous 
satellite the time autocorrelation of the fluctuating signal is given by 


P A (r) ~ exp 



(28) 


Using the same arguments, the cross correlation of signals received at Qj and 
Q will be given by 


where 


P A A (?) exp 
A 1 A 2 



i— ' D/ 


(29) 



Then two antennas spaced £ permit us to determine V , the EW component of 
the drift velocity, by measuring r m . Knowing V D it is then possible to determine 
from the time correlation function. 
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Table 1 

Magnification Factors and Relative Power for 
Satellites at Different Heights 


R i + R s 
Km 

D x 

Km 

M 

Km 

A 2 + R 

Ri s 

A 

36,000 

8 

1 

0 db 

600 

4 

2.2 

-0.8 db 

500 

3.3 

2.5 

-0.9 db 

400 

2.1 

4.0 

-2 db 

325 

0.6 

13 

-20 db 


Note: Quantities in Table I are for the following assumed parameters 


R = 300 km 
= 225 km 
W = 50 km 
f = 136 MHz 

The time r Q = -C x /u D determines the rate of fading of amplitude fluctuations. 

For the case of a satellite at an altitude of the order of 600 km, D x % 4. If a 
circular orbit is assumed, its velocity will be u s % 7.8 x 10 3 m/sec. Since the 
drifting velocity of the ionosphere is in the order of 10 2 m/sec it can be assumed 
that during a time r of the order of (£ x /u s ) the ionosphere is static and the time 
variation of the signals are only due to the satellite movement. As it was pointed 
out this case can be studied by assuming that the transmitter is on the ground 
and that the receiver is at the satellite, then using equation (25) 


where 


p A ( t ) % exp 



(30) 


If equations (19), (25) and (27) are used for a different satellite altitudes Table I 
is obtained. 

This table shows that as the satellite gets near to the layer of irregularities the 
fluctuations of amplitude decrease, but the spacecraft has to be very close to the 
layer for a substantial reduction. 
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DIVERSITY COMBINING TECHNIQUES 


Space diversity techniques can be used to improve the signal to noise ratio 
which suffers degradation by ionospheric irregularities. Diversity combining 
techniques can be studied in a general way to include polarization, frequency, 
space and other types of diversity systems. Here a brief review of this subject 
is made. 

Diversity combining techniques has been studied for many years. Brennan 
(1959) made an analysis of three of the most used techniques: selection diversity, 
maximal ratio diversity and equal gain diversity. 

In selection diversity the best of N noisy signals is chosen and used at a given 
instant of time. 

In maximal ratio diversity the output is given by 


f(0=£\ f. (t) 

i = 1 


where 

f is one of N noisy signals 

q _ rms (signal) 

1 mean square (noise) 

In equal gain diversity we have 


f (t) 


IN 

L 


f. co 


One conclusion from Brennan’s study is that for N = 2 all three types of 
diversity techniques give similar improvement. 


In communication systems the reliability of a link is defined as the percentage 
of time a prescribed signal level is exceeded. Brennan defines diversity gain 
as the reduction in signal level obtained for a diversity system for a given 
reliability. 
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High reliability systems are designed for 99 to 99.9% reliability. If the power 
fluctuations are Rayleigh distributed then the diversity gains shown in Table II 
are obtained. 


Table II 


Diversity Gain for 99% Reliability (Brennan, 1959) 


N 

Selection 

Equal Gain 

Maximal Ratio 

2 

10 

11 

12 

3 

14 

16 

16.5 

4 

16 

18.5 

19.5 


Table II shows that for N = 2 satellite transmitter power reduction of 10 db can 
be obtained for a high reliability link. For lower reliability the diversity gain 
is reduced; at 50% reliability a reduction of power of the order of 3 db is 
possible. 

Table II results were obtained under the assumption that the signals are un- 
correlated and that there is no phase difference between the signals R . 

Brennan found that cross correlations above 0.3 produce negligible results and 
that phase differences below 37.5° produce signal degradations effects. In 
the case under study the antennas should be separated so that the cross corre- 
lations of signals received simultaneously by the two antennas is smaller than 
0.3. Then any of the combining techniques mentioned can be used to improve 
the signal to noise ratio. 

For large separations of the antennas, the problem of phase difference between 
the signals exists. This difference depends on satellite position as shown in 
Figure 3 and is given d sin 0/k if where A-i is the minimum length of the infor- 
mation, is defined by the information bandwidth. For example, if the band- 
width is 300 kHz, d = 1 km then 6 should be greater than 45° to keep the phase 
difference below 37.5°. 

To record the signals from a satellite its angular position has to be determined 
either by manual or automatic track, therefore, the phase difference between 
the signals that arrive at separated antennas can always be approximated. Then, 
degradation caused by phase differences can be reduced to negligible values. 


15 




s 



Figure 3. Signals from infinity received by two spaced antennas have a phase difference 

d sin 8 
Ai 


ANALYSIS OF DATA 

Scintillation data was recorded at Ancon from the beginning of 1970. The data 
was collected for two users, NASA and the Instituto Geofisico del Peru. To date 
there is a total of 45 recordings for the purpose of this experiment. Most of the 
data collected during the months of March and April correspond to one of the 
peaks of scintillation activity (Tisnado, et al, 1970). 

Figure 4 shows a block diagram of the equipment used in this experiment. The 
antennas were oriented in an East- West direction and spaced 366 m. ATS-3 satellite 
transmissions distorted by the ionosphere were received and recorded in mag- 
netic tape. 

The magnetic tapes are then taken to Jicamarca Radar Observatory where the 
computing facilities are used. Figure 5 illustrates the processing system. The 
signals corresponding to the East and West antennas are played back in a FR-600 
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Figure 4. Block diagram of equipment used to gather scintillation data. 
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Figure 5. Block 


of processing equipment for scintillation data. 
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tape recorder and fed to a pair of tracking filters where the signals are syn- 
chronous detected* This type of detection was chosen for its improved sensitivity 
and linear response which permits the observations of low amplitude fluctuations 
with sufficient resolution* The detected signals are then digitized and fed to the 
Jicamarca digital correlator which computes the autocorrelation function corre- 
sponding to each antenna and the cross correlation function between the antennas. 
The resulting information, after 3 minutes of integration time, is fed to the PB- 
250 for normalization and plotting. While the PB-250 is plotting^the computation 
of another set of functions is being done by the correlator, thus one set of func- 
tions can be plotted every 3 minutes. 

The following standard expressions were used to compute the auto- and cross- 
correlation functions 

N 


£ < E , - E > ( E i,i- E ) 

i = 1 



i = 1 1=1 


where 

E T and E n is the signal amplitude of an antenna 
E is the average signal amplitude 

N 

£ (e; - e ') ( E r tJ - I") 

Pa' a" ( r i) =~7 ~r z = 

/N N 

]/£ ( E 'i - e > 2 £ ( E r - e ") 2 

i = 1 i = l 

where E* and E" are signals corresponding to each of the antennas. 

Figure 6 shows three typical correlation functions obtained. They are repre- 
sentative of what has been observed. The fact that the peak of the cross corre- 
lation functions is nearly unity, indicates that the ionospheric irregularities can 
be considered "frozen” as they pass over the antennas. 

The width of the autocorrelation function, Tq , determines the fading rate and the 
scale sizes of the irregularities. The interval r m , at which the maximum of the 
cross correlation occurs, determines the East-West component of the drift 
velocity of irregularities. To obtain this information, Guassian functions were 
fitted to the auto- and cross -correlation functions, then equations (28) and (29) 
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A - Autocorrelation 
B - Crosscorrelation 




Figure 6. Three typical auto and crosscorrelation signals. 




were used. First, V_, the East-West component of the drift velocity is com- 
puted from r m taking V D = 366 rn/r m ; next the correlation length l x is completed 
from t q and V D tracking l x = V D ’ T 0 • 

The result of these computations are presented in Figures 7, 8, and 9 which 
show the distribution of correlation length. West to East drift velocity, and 
correlation width. 

As it was mentioned earlier in this report the height of the irregularities is an 
important parameter in the design of space diversity systems, since the corre- 
lation distance on the ground will depend on it for low orbiting satellites. To 
obtain the height information of the irregularities observations are being con- 
ducted with the Jicamarca Radar. The layered structure of the irregularities 
is clearly established with these observations. Figures 10 and 11 show two 
typical film records obtained with the Jicamarca radar. 

The case of multiple layers is not uncommon but at this time is has not been 
established whether or not all the layers are affecting sate Rite transmissions. 

In many cases Ancon has reported no scintillations while in Jicamarca’s records, 
irregularities were clearly present. 
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Figure 10. Typical radar observations of a layer of irregularities. 
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JICAMA RCA 24 March 1970 



Figure 11. Typical radar observations of multiple layer of irregularities. 




At this time in Jicamarca a relative measurement of power reflected from the 
irregularities is being performed. It is hoped that this will permit a determina- 
tion of the threshold at which satellite transmissions are affected. 

There are two factors that have to be considered when comparing Jicamarca and 
Ancon measurements. First, both stations are not observing the same region of 
the ionosphere. Second, Jicamarca and Ancon are sensitive to different scale 
sizes of irregularities due to different operating frequencies. 


RECOMMENDATIONS FOR A SPACE DIVERSITY SYSTEM 

Although it is still necessary to collect more data to improve the statistics of 
our results, some general recommendation can be given at this time about a 
space diversity system. 

For synchronous satellites the spacing between antennas should be greater than 
300 m to obtain low correlation. The low value of the cross correlation co- 
efficient for t = 0 in Figure 6 indicates that the separation of 366 m now in use 
at Ancon will be adequate for a diversity system. 

For lower satellites the effect of magnification of the sizes determines that 
greater separation should be used. Since most of the satellites have altitudes 
above 500 km the estimated magnification factors are of the order of 3 as ob- 
tained from Table I. This gives separations of the order of 1 km which is not 
at all impractical. Except for practical limitations any reasonable distance 
above 1 km can be used for lower satellites. How far above one kilometer 
would depend on the percentage of reliability expected since there are times 
at which the irregularities exist at all heights in the range between 300 and 
1000 km (Farley, et al, 1970). 

Two practical problems associated with large separations are cable attenuation 
and phase differences between the signals. The problem of cable attenuation can 
be easily solved with the use of preamplifiers at both antennas to compensate for 
the losses. As was explained in this report phase differences due to ionospheric 
fading of satellite signals, in a space diversity system are a function of satellite 
time delay or position. The fact that this position is known permits one to com- 
pensate adequately for the phase differences with reasonable spacing. 

It has been shown (Brennan, 1959) that the phase differences do not need to be com- 
pensated with high precision to permit us to solve this problem. As the satellite 
moves past the station the position information can be fed to a programmed device 
which will introduce in steps the appropriate phase shift or delay. Such a device 
has already been demonstrated by Goddard Space Flight Center. 
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FUTURE WORK 

The height of the irregularities that cause satellite scintillations is an important 
parameter since scale sizes on the ground depend on them. Simultaneous obser- 
vations between Jicamarca and Ancon will be analyzed to obtain this infoi mation. 


Lower satellites will be observed to measure the scale sizes and thus substanti- 
ate our estimations for the magnification parameter. These observations will be 
also used to estimate the height of the irregularities. 

More tapes are being processed and more observations will be conducted during 
the September to December period of high scintillation activity. Time variation 
of drift velocity and scintillation index will be obtained. This will provide im- 
portant information about the phenomenon. 

Propagation through multiple layers will be studied and the distribution of 
amplitude fluctuations will be determined. 
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